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(54) Abstract Title 

Modelling the dynamic behaviour of a complete drilling tool assembly 

(57) The dynamic behaviour of at least a drill bit and a section of drill string are determined using a numerical 
model. Interactions between the drill bit, drill string, rock formation, drilling fluid and the above surface 
assembly are all incorporated, allowing the method to be used to design and optimise drilling systems via an 
iterative approach. The model output may be converted into a visual representation of the drilling process. 
Parameters that can be simulated and used to measure drilling performance include the rate of penetration, 
weight on bit rate of wear of bit and deviation from the desired drill path. The numerical calculations may be 
performed using a finite element model. 
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SIMULATING THE DYNAMIC RESPONSE OF A DRILLING TOOL 
ASSEMBLY AND ITS APPLICATION TO DRILLING TOOL 
ASSEMBLY DESIGN OPTIMIZATION AND DRILLING 
PERFORMANCE OPTIMIZATION 

FTELD OF THE INVENTION 

The invention relates generally to drilling a wellbore, and more specifically 
to simulating the drilling performance of a drilling tool assembly drilling a 
wellbore. In particular, the invention relates to methods for simulating the 
5 dynamic response of a drilling tool assembly, methods for optimizing a drilling 
tool assembly design, and methods for optimizing the drillmg performance of a 
drilling tool assembly. 

BACKGROUND OF THE INVENTION 

Figure 1 shows one example of a conventional drilling system for drilling 
10 an earth formation. The drilling system includes a drilling rig 10 used to turn a 
drilling tool assembly 12 which extends downward into a wellbore 14. The 
drilling tool assembly 12 includes a drilling string 16, and a bottomhole assembly 
(BHA) 18, attached to the distal end of the drill string 16. 

The drill string 16 comprises several joints of drill pipe 16a connected end 
15 to end through tool joints 16b. The drill string 16 transmits drilling fluid (through 
its hollow core) and transmits rotational power from the drill rig 10 to the BHA 
18. In some cases the drill string 16 further includes additional components such 
as subs, pup joints, etc. 
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The BHA 18 includes at least a drill bit 20. Typical BHAs may also 
include additional components attached between the drill string 16 and the drill bit 
20. Examples of additional BHA components include drill collars, stabilizers, 
measurement-while-drilling (MWD) tools, logging-while-drilling (LWD) tools, 
5 and downhole motors. 

In general, drilling tool assemblies 12 may include other drilling 
components and accessories, such as special valves, such as kelly cocks, blowout 
preventers, and safety valves. Additional components included in a drilling tool 
assembly 12 may be considered a part of the drill string 16 or a part of the BHA 18 
1 0 depending on their locations in the drilling tool assembly 12. 

The drill bit 20 in the BHA 18 may be any type of drill bit suitable for 
drilling earth formation. Two common types of earth boring bits used for drilling 
earth formations are fixed-cutter (or fixed-head) bits and roller cone bits. Figure 2 
shows one example of a fixed-cutter bit. Figure 3 shows one example of a roller 
15 cone bit. 

Referring to Figure 2, fixed-cutter bits (also called drag bits) 21 typically 
comprise a bit body 22 having a threaded connection at one end 24 and a cutting 
head 26 formed at the other end. The head 26 of the fixed-cutter bit 21 typically 
comprises a plurality of ribs or blades 28 arranged about the rotational axis of the 

20 bit and extending radially outward firom the bit body 22. Cutting elements 29 are 
embedded in the raised ribs 28 to cut formation as the bit is rotated on a bottom 
surface of a wellbore. Cutting elements 29 of fixed-cutter bits typically comprise 
polycrystalline diamond compacts (PDC) or specially manufactured diamond 
cutters. These bits are also referred to as PDC bits. 

25 Referring to Figure 3, roller cone bits 30 typically comprise a bit body 32 

having a threaded connection at one end 34 and a plurality of legs (not shown) 
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extending from the other end. A roller cone 36 is mounted on each of the legs and 
is able to rotate with respect to the bit body 32. On each cone 36 of the bit 30 are 
a plurality of cutting elements 38, typically airanged in rows about the surface of 
the cone 36 to contact and cut through formation encountered by the bit. Roller 

5 cone bits 30 are designed such that as a drill bit rotates, the cones 36 of the bit 30 
roll on the bottom surface of the wellbore (called the "bottomhole") and the 
cutting elements 38 scrape and crush the formation beneath them. In some cases, 
the cutting elements 38 on the roller cone bit 30 comprise milled steel teeth 
formed on the surface of the cones 36. In other cases, the cutting elements 38 

10 comprise inserts embedded in the cones. Typically, these inserts are tungsten 
carbide inserts or polycrystalline diamond compacts. In some cases hardfacing is 
applied to the surface of the cutting elements to improve wear resistance of the 
cutting structure. 

For a drill bit 20 to drill through formation, sufficient rotational moment 
15 and axial force must be applied to the bit 20 to cause the cutting elements of the 
bit 20 to cut into and/or crush formation as the bit is rotated. The axial force 
applied on the bit 20 is typically referred to as the **weight on bit" (WOB). The 
rotational moment applied to the drilling tool assembly 12 at the drill rig 10 
(usually by a rotary table) to turn the drilling tool assembly 12 is referred ^o as the 
20 "rotary torque". The speed at which the rotary table rotates the drilling tool 
assembly 12, typically measured in revolutions per minute (RPM), is referred to as 
the "rotary speed". Additionally, the portion of the weight of the drilling tool 
assembly supported at the rig 10 by the suspending mechanism (or hook) is 
typically referred to as the hook load. 
25 During drilling, the actual WOB is not constant. Some of the fluctuation in 

the force applied to the bit may be the result of the bit contacting with formation 
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having harder and softer portions that break unevenly. However, in most cases, 
the majority of the fluctuation in the WOB can be attributed to driUing tool 
assembly vibrations. Drilling tool assemblies can extend more than a mile in 
length while being less than a foot in diameter. As a result, these assemblies are 
5 relatively flexible along their length and may vibrate when driven rotationally by 
the rotary table. Several modes of vibration are possible for .drilling tool 
assemblies. In general, drilling tool assemblies may experience torsional, axial 
and lateral vibrations. Although partial damping of vibration may result due to 
viscosity of drilling fluid, friction of the drill pipe robbing against the wall of the 

10 wellbore, energy absorbed in drilling the formation, and drilling tool assembly 
impacting with wellbore wall, these sources of damping are typically not enough 
to suppress vibrations completely. 

Up to now, vibrations of a drilling tool assembly have been difficult to 
predict because different forces may combine to produce the various modes of 

15 vibration, and models for simulating the response of an entire drilling tool 
assembly including roller cone bit interacting with formation in a drilling 
environment have not been available. However, drilling tool assembly vibrations 
are generally undesirable, not only because they are difficult to predict, but also 
because they can significantly affect the instantaneous force applied on the bit. 

20 This can result in the bit not operating as expected. For example, vibrations can 
result in off-centered drilling, slower rates of penetration, excessive wear of the 
cutting elements, or premature failure of the cutting elements and the bit. Lateral 
vibration of the drilling tool assembly may be a result of radial force imbalances, 
mass imbalance, and bit/formation interaction, among other things. Lateral 

25 vibration results in poor drilling tool assembly performance, overgage hole 
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drilling, out-of-round, or "lobed" wellbores and premature failure of both the 
cutting elements and bit bearings. 

When the bit wears out or breaks during drilling, the entire drilling tool 
assembly must be lifted out of the wellbore section-by-section and disassembled in 
5 an operation called a "pipe trip". In this operation, a heavy hoist is required to pull 
the drilling tool assembly out of the wellbore in stages so that each stand of pipe 
(typically pipe sections of about 90 feet) can be unscrewed and racked for the later 
re-assembly. Because the length of a drilling tool assembly may extend for more 
than a mile, pipe trips can take several hours and can pose a significant expense to 

10 the wellbore operator and drilling budget. Therefore, the ability to design drilling 
tool assemblies which have increased durability and longevity, for example, by 
minimizing the wear on the drilling tool assembly due to vibrations, is very 
important and greatly desired to minimize pipe trips out of the wellbore and to 
more accurately predict the resulting geometry of the wellbore drilled. 

15 Simulation methods have been previously introduced which characterize 

either the interaction of a bit with the bottomhole surface of a wellbore or the 
dynamics of a bottomhole assembly (BHA). However, no prior art simulation 
techniques have been developed to cover the dynamic modeling of an entire 
drilling tool assembly. As a result, the dynamic response of a drilling tool 

20 assembly or the effect of a change in configuration on drilling tool assembly 
performance can not be accurately predicted. 

One simulation method for characterizing interaction between a roller cone 
bit and an earth formation is described in U.S. Patent Application No. 09/524,088, 
entitled "Method for Simulating Drilling of Roller Cone Bits and its Application to 

25 Roller Cone Bit Design and Perforaiance", and assigned to the assignee of the 
present invention. This application discusses general methods for predicting 
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cutting element interaction with earth formations. The application also discussed 
types of experimental tests that can be performed to obtain cutting 
element/formation interaction data. Another simulation method for characterizing 
cutting element/fomiation interaction for a roller cone bit is described in Society 
5 of Petroleum Engineers (SPE) Paper No. 29922 by D. Ma et ah, entitled, "The 
Computer Simulation of the Interaction Between Roller Bit and Rock?. 

Methods for optimizing tooth orientation on a roller cone bits are disclosed 
in PCT International Publication No. WOOO/12859 entitled, "Force-Balanced 
Roller-Cone Bits, Systems, Drilling Methods, and Design Methods*' and PCT 

10 International Publication No. WOOO/12860 entitled, "Roller-Cone Bits, Systems, 
Drilling Methods, and Design Metfiods with Optimization of Tooth Orientation. 

Similarly, SPE Paper No. 15618 by T. M. Warren et aL, entitled "Drag Bit 
Performance Modeling*' discloses a method for simulating the performance of 
PDC bits. Also disclosed are methods for defining the bit geometry, and methods 

IS for modeling forces on cutting elements and cutting element wear during drilling 
based on experimental test data. Examples of experimental tests that can be 
performed to obtain cutting element/earth formation interaction data are also 
disclosed. Experimental methods that can be performed on bits in earth 
formations to characterize bit/earth formation interaction are discussed in SPE 

20 Paper No. 15617 by T. M. Warren et al., entitled "Laboratory Drilling 
Performance of PDC Bits". 

While prior art simulation methods, such as those described above cover 
either the interaction of the bit with the formation or the BHA dynamics, no prior 
art simulation technique has been developed to cover the dynamic modeling of the 

25 entire drilling tool assembly. As a result, accurately predicting the response of a 
drilling tool assembly has been virtually impossible. Additionally, the change in 
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the dynamic response of a drilling tool assembly when a component of the drilling 
tool assembly is changed is not well understood. 

In view of the above it is clear that a method for simulating the dynamic 
response of an entire drilling tool assembly, which takes into account bit 
5 interaction with the bottom surface of the wellbore, drilling tool assembly 
interaction with the wall of the wellbore and damping effects of the drilling fluid 
on the drill pipe, is both needed and desired. Additionally, a model for predicting 
changes in drilling tool assembly performance due to changes in drilling tool 
assembly , configuration, and for determining optimal drilling tool assembly 
10 designs and/or optimal drilling operating parameters (WOB, RPM, etc.) for a 
particular depth, fomiation, and/or drilling tool assembly is desired. 

SUMMARY OF THE INVENTJON 

The invention provides methods for simulating the dynamic response of a 
drilling tool assembly drilling an earth formation. The drilling tool assembly 

1 5 comprises at least a drill pipe and a drill bit. Methods for simulating the dynamic 
response of drilling tool assemblies may be used to generate a visual 
representation of drilling, to design drilling tool assemblies, and to optimize the 
drilling performance of a drilling tool assembly. 

One method for generating a visual representation of a drilling tool 

20 assembly which comprising at least a drill pipe and a drill bit comprises solving 
for a dynamic response of the drilling tool assembly to an incremental rotation, 
determining, based on the dynamic response, parameters of craters removed from 
a bottomhole surface of the formation due to contact of the bit with the bottomhole 
surface during the incremental rotation, and calculating a bottomhole geometry, 

25 wherein the craters are removed from the bottomhole surface. The method further 
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comprises repeating the solving, determining, and calculating for a select number 
of successive incremental rotations, and converting the dynamic responses and the 
bottomhole geometry parameters into a visual representation. 

One method for optimizing a drilling tool assembly design comprises 
5 simulating a dynamic response of the drilling tool assembly, adjusting a value of at 
least one drilling tool assembly design parameter, and repeating the simulating. 
The method further comprises repeating the adjusting and the simulating until at 
least one drilling performance parameter is determined to be at an optimum value. 
One method for determining at least one optimal drilling operating 

10 parameter for a drilling tool assembly comprises simulating a dynamic response of 
the drilling tool assembly, adjusting the value of at least one drilling operating 
parameter, and repeating the simulating. The method further includes repeating 
the adjusting and the simulating until at least one drilling performance parameter 
is determined to be at an optimal value. 

15 One method for designing a drilling tool assembly comprises defining 

initial drilling tool assembly design parameters, simulating the dynamic response 
of the drilling tool assembly, adjusting a value of at least one of the drilling tool 
assembly design parameters, and repeating the simulating and the adjusting a 
select number of times. The method further comprises evaluating the dynamic 

20 responses, and selecting, based on the evaluating, desired drilling tool assembly 
design parameters. 

Other aspects and advantages of the invention will be apparent from the 
following description and the appended claims. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



Figure 1 shows a schematic diagram of a prior art drilling system for 
drilling earth formations. 

Figure 2 shows a perspective view of a prior art fixed-cutter bit. 
5 Figure 3 shows a perspective view of a prior art roller cone bit 

Figure 4, shows one example of drilling tool assembly. 

Figure 5 shows a flow chart of one embodiment of a method for simulating 
the dynamic response of a drilling tool assembly. 

Figure 6 shows a flow chart of one method of incrementally solving for the 
1 0 dynamic response of a drilling tool assembly. 

Figures 7A-C shows a more detailed flow chart of a method for 
incrementally solving for the dynamic response of a drilling tool assembly in 
which constraint loads are updated to account for interaction between the drilling 
tool assembly and the drilling environment during the incremental rotation. 
1 5 Figure 8 shows a general flow chart of one method for determining an 

optimal value of at least one drilling tool assembly design parameter. 

Figure 9 shows a more detailed flow chart of a method for determining an 
optimal value of at least one drilling tool assembly design parameter. 

Figure 10 shows a general flow chart of one method for determining an 
20 optimal value for at least one drilling operating parameter for a drilling tool 
assembly. 

Figure 1 1 shows a more detailed flow chart of a method for determining an 
optimal value for at least one drilling operating parameter for a drilling tool 
assembly. 

25 Figure 12 shows one example of output data converted into a visual 

representation. 
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DETAILED DFSCRTPTION OF THF INVENTION 



The invention provides methods for simulating the dynamic response of a 
drilling tool assembly drilling an earth formation, methods for optimizing a 
drilling tool assembly design, and methods for optimizmg drilling tool assembly 
5 performance. 

In accordance with the invention, a drilling tool assembly comprises at least 
one segment (or joint) of drill pipe and a drill bit. The components of a drilling 
tool assembly may be more generally referred to as a drill string and a bottomhole 
assembly (BHA). The drill string comprises one or more joints of drill pipe. The 

1 0 BHA comprises at least a drill bit. 

In a typical drilling tool assembly, the drill string comprises several joints 
of drill pipe connected end to end, and the bottomhole assembly comprises one or 
more drill collars and a drill bit attached to an end of the BHA. The drill string 
may further include additional components, such as a kelly, kelly cocks, blowout 

15 preventers, safety valves, etc. The BHA may forther include additional 
components, such as stabilizers, a downhole motor, MWD tools, and LWD tools, 
for example. Therefore, in accordance with the invention, a drilling tool assembly 
may be as simple as a single segment of drill pipe attached to a drill bit, or as 
complex as a multi-component drill string which includes a kelly, a lower kelly 

20 cock, a kelly saver sub, several joints of drill pipe with tool joints, etc., and a 
multi-component BHA which includes drill collars, stabilizers, and additional 
specialty items (e.g., subs, pup joints, reamers, valves, MWD tools, LWD tools, 
and a drill bit). 

While in practice, a BHA comprises at least a drill bit, in embodiments of 
25 the invention discussed below, the parameters of the drill bit, required for 
modeling interaction between the drill bit and the bottomhole surface, are 



generally considered separately firom the BHA parameters. This separate 
consideration of the bit allows for interchangeable use of any drill bit model as 
determined by the system designer. 

One example of a drilling tool assembly 50 is shown in Figure 4. In this 
5 embodiment, the drilling tool assembly is suspended from a hook 62 and rotated 
by a rotary table 64. The drilling tool assembly 50 comprises a drill string 52 and 
BHA 54. The drill string 52 comprises a plurality of joints of drill pipe 56. The 
BHA 54 comprises a drill collar 58 and a drill bit 60. The drill bit 62 shown in 
this example is a roller cone drill bit. In other embodiments any type of drill bit 
10 may be used. 

To simulate the dynamic response of a drilling tool assembly, such as the 
one shown in Figure 4, for example, components of the drilling tool assembly 
need to be mathematically defined. For example, the drill string may generally be 
defined in terms of geometric and material parameters, such as the total length, the 

15 total weight, inside diameter (ID), outside diameter (OD), and material properties 
of the various components of the drill string. Material properties of the drill string 
components may include the strength, and elasticity of the component material. 
Each component of the drill string may be individually defined or various parts 
may be defined in the aggregate. For example, a drill string comprising a plurality 

20 of substantially identical joints of drill pipe may be defined by the number of drill 
pipe joints of the drill string, and the ID, OD, length, and material properties for 
one drill pipe joint. Similarly, the BHA may be defined in terms of parameters, 
such as the ID, OD, length, and material properties of one drill collar and of any 
other component that makes up the BHA. 

25 The geometry and material properties of the drill bit also need to be defined 

as required for the method selected for simulating drill bit interaction with the 
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earth fonnation at the bottom surface of the wellbore. One example of a method 
for simulating a roller cone drill bit drilling an earth formation can be found in the 
previously mentioned U.S. Patent Application No. 09/524,088, assigned to the 
assignee of the present invention and now incorporated herein by reference in its 
5 entirety. 

To simulate the dynamic response of a drilling tool assembly drilling earth 
formation, the wellbore trajectory, in which the drilling tool assembly is to be 
confined also needs to be defined along with an initial wellbore bottom surface 
geometry. Because the wellbore trajectory may be straight, curved, or a 

10 combination of straight and curved sections, wellbore trajectories, in general, may 
be defined by defining parameters for each segment of the trajectory. For 
example, a wellbore comprising N segments may be defined by the length, 
diameter, incHnation angle, and azimuth direction of each segment and an 
indication of the order of the segments (/.e., first, second, etc.). Wellbore 

15 parameters defined in this manner can then be used to mathematically produce a 
model of the entire wellbore trajectory. Fonnation material properties along the 
wellbore may also be defined and used. Additionally, drilling operating 
parameters, such as the speed at which the drilling tool assembly is rotated and the 
hook load (weight of the drilling tool assembly suspended at the hook 62), also 

20 need to be defined. 

Once the parameters of the system (drilling tool assembly under drilling 
conditions) are defined, they can be used along with various interaction models to 
simulate the dynamic response of the drilling tool assembly drilling earth 
formation as described below. 

25 

Method for Simulating the Dynamic Response of Drilling Tool Assemblv 
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In one aspect, the invention provides a method for simulating the dynamic 
response of a drilling tool assembly drilling earth formation. Advantageously, this 
method takes into account interaction between the entire drilling tool assembly 
and the drilling environment. Interaction between the drilling tool, assembly and 
5 the drilling environment may include interaction between the drill bit at the end of 
the drilling tool assembly and the formation at the bottom of the wellbore. 
Interaction between the drilling tool assembly and the drilling environment also 
may include interaction between the drilling tool assembly and the side (or wall) 
of the wellbore. Further, interaction between the drilling tool assembly and 

10 drilling environment may include viscous damping effects of the drilling fluid on 
the dynamic response of the drilling tool assembly. 

A flow chart for one embodiment of the invention is illustrated in Figure 5. 
The first step in this embodiment is selecting (defining or otherwise providing) 
parameters 100, including initial drilling tool assembly parameters 102, initial 

15 drilling environment parameters 104, drilling operating parameters 106, and 
drilling tool assembly/drilling environment interaction information (parameters 
and/or models) 108. The nest step involves constructing a mechanics analysis 
model of the drilling tool assembly 110. The mechanics analysis model can be 
constructed using the drilling tool assembly parameters 102 and Newton's law of 

20 motion. The next step involves determining an initial static state of the drilling 
tool assembly 112 in the selected drilling environment using the mechanics 
analysis model 110 along with drilling environment parameters 104 and drilling 
tool assembly/drilling environment interaction information 108. Once the 
mechanics analysis model is constructed and an initial static state of the drill string 

25 is determined, the resulting static state parameters can be used with the drilling 
operating parameters 106 to incrementally solve for the dynamic response 114 of 
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the drilling tool assembly 50 to rotational input from the rotary table 64 and the 
hook load provided at the hook 62. Once a simulated response for an increment in 
time (or for the total time) is obtained, results from the simulation can be provided 
as output 118, and used to generate a visual representation of drilling if desired. 
5 In one example, illustrated in Figure 6, incrementally solving for the 

dynamic response (indicated as 116) may not only include solving the mechanics 
analysis model ifor the dynamic response to an incremental rotation, at 120, but 
may also include determining, from the response obtained, loads (e.g., drilling 
environment interaction forces) on the drilling tool assembly due to interaction 

10 between the drilling tool assembly and the drilling environment during the 
incremental rotation, at 122, and resolving for the response of the drilling tool 
assembly to the incremental rotation, at 124, under the newly determined loads. 
The determining and resolving may be repeated in a constraint update loop 128 
until a response convergence criterion 126 is satisfied. Once a convergence 

15 criterion is satisfied, the entire incremental solving process 116 may be repeated 
for successive increments until an end condition for simulation is reached. 

During the simulation, the constraint forces initially used for each new 
incremental calculation step may be the constraint forces determined during the 
last incremental rotation. In the simulation, incremental rotation calculations are 

20 repeated for a select number of successive incremental rotations until an end 
condition for simulation is reached. A more detailed example of an embodiment 
of the invention is shown in Figure 7A-C. 

For the example shown in Figure 7A-C, the parameters provided as input 
200 include drilling tool assembly design parameters 202, initial drilling 

25 environment parameters 204, drilling operating parameters 206, and drilling tool 
assembly/drilling environment interaction parameters and/or models 208. 
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Drilling tool assembly design parameters 202 may include drill string 
design parameters, BHA design parameters, and drill bit design parameters. In the 
example shown, the drill string comprises a plurality of joints of drill pipe, and the 
BHA comprises drill collars, stabilizers, bent housings, and other downhole tools 
5 (e.g., MWD tools, LWD tools, downhole motor, etc.), and a drill bit. As noted 
above, while the drill bit, generally, is considered a part of the BHA, in this 
example the design parameters of the drill bit are shown separately to illustrate 
that any type of drill bit may be defined and modeled using any drill bit analysis 
model. 

10 Drill string design parameters include, for example, the length, inside 

diameter (ID), outside diameter (OD), weight (or density), and other material 
properties of tibe drill string in the aggregate. Alternatively, drill string design 
parameters may include the properties of each component of the drill string and 
the number of components and location of each component of the drill string. For 

1 5 example, the length, ID, OD, weight, and material properties of one joint of drill 
pipe may be provided along with the number of joints of drill pipe which make up 
the drill string. Material properties used may include the type of material and/or 
the strength, elasticity, and density of the material. The weight of the drill string, 
or individual components of the drill string may be provided as "weight in drilling 

20 fluids" (the weight of the component when submerged in the selected drilling 
fluid). 

BHA design parameters include, for example, the bent angle arid 
orientation of the motor, the length, equivalent inside diameter (ID), outside 
diameter (OD), weight (or density), and other material properties of each of the 
25 various components of the BHA. In this example, the drill collars, stabilizers, and 
other downhole tools are defined by their lengths, equivalent IDs, ODs, material 
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properties, weight in drilling fluids, and position in the drilling tool assembly. 

The drill bit design parameters include, for example, the bit type (roller 
cone, fixed-cutter, etc.) and geometric parameters of the bit. Geometric 
parameters of the bit may include the bit size (e.g., diameter), number of cutting 

5 elements, and the location, shape, size, and orientation of the cutting elements. In 
the case of a roller cone bit, drill bit design parameters may further include cone 
profiles, cone axis offset (offset fi-om perpendicular with the bit axis of rotation), 
the number of cutting elements on each cone, the location, size, shape, orientation, 
etc. of each cutting element on each cone, and any other bit geometric parameters 

1 0 {e.g., journal angles, element spacings, etc.) to completely define the bit geometry. 
In general, bit, cutting element, and cone geometry may be converted to 
coordinates and provided as input. One preferred method for obtaining bit design 
parameters is the use of 3-dimensional CAD solid or surface models to facilitate 
geometric input. Drill bit design parameters may further include material 

1 5 properties, such as strength, hardness, etc. of components of the bit. 

Initial drilling environment parameters 204 include, for example, wellbore 
parameters. Wellbore parameters may include wellbore trajectory (or geometric) 
parameters and wellbore formation parameters. Wellbore trajectory parameters 
may include an initial wellbore measured depth (or length), wellbore diameter, 

20 inclination angle, and azimuth direction of the wellbore trajectory. In the typical 
case of a wellbore comprising segments having different diameters or differing in 
direction, the wellbore trajectory information may include depths, diameters, 
inclination angles, and azimuth directions for each of the various segments. 
Wellbore trajectory information may further include an indication of the curvature 

25 of the segments (which may be used to determine the order of mathematical 
equations used to represent each segment). Wellbore formation parameters may 
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include the type of formation being drilled and/or material properties of the 
formation such as the formation strength, hardness, plasticity, and elastic modulus. 

Drilling operating parameters 206, in this embodiment, include the rotary 
table speed at which the drilling tool assembly is rotated (RPM), the downhole 
5 motor speed if a downhole motor is included, and the hook load. Drilling 
operating parameters 206 may further include drilling fluid parameters, such as the 
viscosity and density of the drilling fluid, for example. It should be understood 
that drilling operating parameters 206 are not limited to these variables. In other 
embodiments, drilling operating parameters 206 may include other variables, such 

10 as, for example, rotary torque and drilling fluid flow rate. Additionally, drilling 
operating parameters 206 for the purpose of simulation may further include the 
total number of bit revolutions to be simulated or the total drilling time desired for 
simulation. However, it should be understood that total revolutions and total 
drilling time are simply end conditions that can be provided as input to control the 

15 stopping point of simulation, and are not necessary for the calculation required for 
simulation. Additionally, in other embodiments, other end conditions may be 
provided, such as total drilling depth to be simulated, or by operator command, for 
example. 

Drilling tool assembly/drilling environment interaction information 208 
20 * includes, for example, cutting element/earth formation interaction models (or 
parameters) and drilling tool assembly/formation impact, friction, and damping 
models and/or parameters. Cutting element/earth formation interaction models 
may include vertical force-penetration relations and/or parameters which 
characterize the relationship between the axial force of a selected cutting element 
25 on a selected formation and the corresponding penetration of the cutting element 
into the formation. Cutting element/earth formation interaction models may also 
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include lateral force-scraping relations and/or parameters which characterize the 
relationship between the lateral force of a selected cutting element on a selected 
formation and the corresponding scraping of the formation by the cutting element. 
Cutting element/formation mteraction models may also include brittle fracture 

5 crater models and/or parameters for predicting formation craters which will likely 
result in brittle fracture, wear models and/or parameters for predicting cutting 
element wear resulting from contact with the formation, and cone shell/formation 
or bit body/formation interaction models and/or parameters for determining forces 
on the bit resulting from cone shell/formation or bit body/formation interaction. 

10 One example of methods for obtaining or determining drilling tool 
assembly/foimation interaction models or parameters can be found in previously 
noted U.S. Patent Application No. 09/524,088, assigned to the assignee of the 
present invention and incorporated herein by reference. Other methods for 
modeling drill bit interaction with a formation can be found in the previously 

15 noted SPE Papers No. 29922, No. 1 5617, and No. 1 561 8, and PCX International 
Publication Nos. WO 00/12859 and WO 00/12860. 

Drilling tool assembly/formation impact, fnction, and damping models 
and/or parameters characterize impact and friction on the drilling tool assembly 
due to contact with the wall of the wellbore and the viscous damping effects of the 

20 drilling fluid. These models/parameters include, for example, drill string- 
BHA/formation impact models and/or parameters, bit body/formation impact 
models and/or parameters, drill string-BHA/formation fhction models and/or 
parameters, and drilling fluid viscous damping models and/or parameters. One 
skilled in the art will appreciate that impact, fnction and damping 

25 models/parameters may be obtained through laboratory experimentation, in a 
method similar to that disclosed in the prior art for drill bits interaction 
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models/parameters. Alternatively, these models may also be derived based on 
mechanical properties of the fomiation and the drilling tool assembly, or may be 
obtained from literature. Prior art methods for determining impact and friction 
models are shown, for example, in papers such as the one by Yu Wang and 
5 Matthew Mason, entitled *Two-Dimensional Rigid-Body Collisions with 
Friction", Journal of Applied Mechanics, Sept. 1992, Vol 59, pp. 635r642. 

As shown in Figure 7A-C, once input parameters/models 200 are selected, 
detemiined, or otherwise provided, a two-part mechanics analysis model of the 
drilling tool assembly is constmcted (at 210) and used to determine the initial 

10 static state (at 232) of the drilling tool assembly in the wellbore. The first part of 
the mechanics analysis model 210a takes into consideration the overall structure 
of the drilling tool assembly, with the drill bit being only generally represented. In 
this embodiment, for example, a finite element method is used (generally 
described at 212) wherein an arbitrary initial state (such as hanging in the vertical 

15 mode free of bending stresses) is defined for the drilling tool assembly as a 
reference and the drilling tool assembly is divided into N elements of specified 
element lengths (z.e., meshed). The static load vector for each element due to 
gravity is calculated. Then element stiffiiess matrices are constructed based on the 
material properties (e.g., elasticity), element length, and cross sectional 

20 geometrical properties of drilling tool assembly components provided as input and 
are used to constroct a stiffiiess matrix, at 212, for the entire drilling tool assembly 
(wherein the drill bit is generally represented by a single node). Similarly, 
element mass matrices are constructed by determining the mass of each element 
(based on material properties, etc.) and are used to construct a mass matrix, at 214, 

25 for the entire drilling tool assembly. Additionally, element damping matrices can 
be constructed (based on experimental data, approximation, or other method) and 
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used to construct a damping matrix, at 216, for the entire drilling tool assembly. 
Methods for dividing a system into finite elements and constructing corresponding 
stif&ess, mass, and damping matrices are known in the art and thus are not 
explained in detail here. Examples of such methods are shown, for example, in 
5 "Finite Elements for Analysis and Design" by J. E. Akin (Academic Press, 1 994). 

The second part 210b of the mechanics analysis model 210 of the drilling 
tool assembly is a mechanics analysis model of tiie drill bit 210b which takes into 
account details of selected drill bit design. The drill bit mechanics analysis model 
210b is constructed by creating a mesh of the cutting elements and cones (for a 

10 roller cone bit) of the bit, and establishing a coordinate relationship (coordinate 
system transformation) between the cutting element and the cones, between the 
cones and the bit, and between the bit and the tip of the BHA. As previously 
noted, examples of methods for constructing mechanics analysis models for roller 
cone drill bits can be found in U.S. Patent Application No. 09/524,088, as well as 

1 5 SPE Paper No. 29922, and PCT International Publication Nos. WO 00/12859 and 
WO 00/12860, noted above. 

Because the response of the drilling tool assembly is subject to the 
constraint within the wellbore, wellbore constraints for the drilling tool assembly 
are determined, at 222, 224. Furst, the trajectory of the wall of the wellbore, which 

20 constrains the drilling tool assembly and forces it to conform to the wellbore path, 
is constructed at 220 using wellbore trajectory parameters provided as input at 
204. For example, a cubic B-spline method or other interpolation method can be 
used to approximate wellbore wall coordinates at depths between the depths 
provided as input data. The wall coordinates are then discretized (or meshed), at 

25 224 and stored. Similarly, an initial wellbore bottom surface geometry, which is 
either selected or determined, is also be discretized, at 222, and stored. The initial 
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bottom surface of the wellbore may be selected as flat or as any other contour, 
which can be provided as wellbore input at 204 or 222. Alternatively, the initial 
bottom surface geometry may be generated or approximated based on the selected 
bit geometry. For example, the initial bottomhole geometry may be selected from 
5 a "library" (i.e., database) containing stored bottomhole geometries resulting from 

the use of various bits. 

In this embodiment, a coordinate mesh size of 1 millimeter is selected for 
the wellbore surfaces (wall and bottomhole); however, the coordinate mesh size is 
not intended to be a limitation on the invention. Once meshed and stored, the 

10 wellbore wall and bottomhole geometry, together, comprise the initial wellbore 
constraints within which the drilling tool assembly must operate, thus, within 
which the drilling tool assembly response must be constrained. 

As shown in Figure 7A-C, once the (two-part) mechanics analysis model 
for the drilling tool assembly is constructed 210 (using Newton's second law) and 

15 the wellbore constraints are specified 222, 224, the mechanics model and 
constraints can be used to determine the constraint forces on the drilling tool 
assembly when forced to the wellbore trajectory and bottomhole from its original 
"stress free" state. In this embodiment, the constraint forces on the drilling tool 
assembly are determined by first displacing and fixing the nodes of the drilling 

20 tool assembly so the centcrline of the drilling tool assembly corresponds to the 
centerline of the wellbore, at 226. Then, the corresponding constraining forces 
required on each node (to fix it in this position) are calculated at 228 from the 
fixed nodal displacements using the drilling tool assembly (i.e., system or global) 
stiffness matrix from 212. Once the "centerline" constraining forces are 

25 determined, the hook load is specified, and initial wellbore wall constraints and 
bottomhole constraints are introduced at 230 along the drilling tool assembly and 
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at the bit (lowest node). The centerline constraints are used as the wellbore wall 
constraints. The hook load and gravitational force vector are used to determine the 
WOB. 

As previously noted, the hook load is the load measured at the hook from 
5 which the drilling tool assembly is suspended. Because the weight of the drilling 
tool assembly is known, the bottomhole constraint force (/.e, WOB) can be 
determined as the weight of the drilling tool assembly minus the hook load and the 
frictional forces and reaction forces of the hole wall on the drilling tool assembly. 
Once the initial loading conditions are introduced, the "centerline" 

10 constraint forces on all of the nodes are removed, a gravitational force vector is 
applied, and the static equilibrium position of the assembly within the wellbore is 
determined by iteratively calculating the static state of the drilling tool assembly 
232. Iterations are necessary since the contact points for each iteration may be 
different. The convergent static equilibrium state is reached and the iteration 

15 process ends when the contact points and, hence, contact forces are substantially 
the same for two successive iterations. Along wifli flie static equilibrium position, 
the contact points, contact forces, friction forces, and static WOB on the drilling 
tool assembly are determined. Once the static state of the system is obtained (at 
232) it can be used as the staring point (initial condition) 234 for simulation of the 

20 dynamic response of the drilling tool assembly drilling earth formation. 

As shown in Figure 7A-C, once input data are provided and the static state 
of the drilling tool assembly in the wellbore is determined, calculations in the 
dynamic response simulation loop 240 can be carried out. Briefly summarizing 
the functions performed in the dynamic response loop 240, the drilling tool 

25 assembly drilling earth formation is simulated by "rotating" the top of the drilling 
tool assembly (and the downhole motor, if used) through an incremental angle (at 
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242), and then calculating the response of tiie drilling tool assembly under the 
previously determined loading conditions 244 to the rotation(s). The constraint 
loads on the drilling tool assembly resulting from interaction with the wellbore 
wall during the incremental rotation are iteratively determined (in loop 245) and 
5 are used to update the drilling tool assembly constraint loads global load 
vector), at 248, and the response is recalculated under the updated loading 
condition- The new response is then rechecked to determine if wall constraint 
loads have changed and, if necessary, wall constraint loads are re-determined, the 
load vector updated, and a new response calculated. Then the bottomhole 

10 constraint loads resulting from bit interaction with the formation during the 
incremental rotation are evaluated based on the new response (loop 252), the load 
vector is updated (at 279), and a new response is calculated (at 280). The wall and 
bottomhole constraint forces are repeatedly updated (in loop 285) until 
convergence of a dynamic response solution is determined (Le,, changes in the 

15 wall constraints and bottomhole constraints for consecutive solutions are 
determined to be negligible). The entire dynamic simulation loop 240 is then 
repeated for successive incremental rotations until an end condition of the 
simulation is reached (at 290) or until simulation is otherwise terminated. A more 
detailed description of the elements in the simulation loop 240 follows. 

20 Prior to the start of the simulation loop 240, drilling operating parameters 

206 are specified. As previously noted, the drilling operating parameters 206 
include the rotary table speed, downhole motor speed (if included in the BHA), 
and the hook load. In this example, the end condition for simulation is also 
provided at 204, as either the total number of revolutions to be simulated or the 

25 total time for the simulation. Additionally, the incremental step desired for 
calculations should be defined, selected, or otherwise provided. In the 
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embodiment shown, an incremental time step of At=10"' seconds is selected. 
However, it should be understood that the incremental time step is not intended to 
be a limitation on the invention. 

Once the static state of the system is known (from 232) and tiie operational 
5 parameters are provided, the dynamic response simulation loop 240 can begin. In 
the first step of the simulation loop 240, the current time increment is calculated at 
241, wherein tj+i = ti + At Then, the incremental rotation which occurs during that 
time increment is calculated, at 242. In this embodiment, the formula used to 
calculate an incremental rotation angle at time tj+i is Oi+i=Oi + RPM *At *60, 

1 0 wherein RPM is the rotational speed (in RPM) of the rotary table provided as input 
data (at 204). The calciilated incremental rotation angle is applied proximal to the 
top of the drilling tool assembly (at the node(s) corresponding to the position of 
the rotary table). If a downhole motor is included in the BHA, the downhole 
motor incremental rotation is also calculated and applied to the corresponding 

15 nodes. 

Once the incremental rotation angle and current time are determined, the 
system's new configuration (nodal positions) under the extant loads and the 
incremental rotation is calculated (at 244) using mechanics analysis model 
modified to include the rotational input as an excitation. For example, a direct 

20 integration scheme can be used to solve the resulting dynamic equilibrium 
equations (modified mechanics analysis model) for the drilling tool assembly. The 
dynamic equilibrium equation (like the mechanics analysis equation) can be 
derived using Newton's second law of motion, wherein the constructed drilling 
tool assembly mass, stiffiiess, and damping matrices along with the calculated 

25 static equilibrium load vector can be used to determine the response to the 
incremental rotation. For the example shown in Figure 7A-C, it should be 



24 



understood that at the first time increment ti the extant loads on the system are the 
static equilibrium loads (calculated for to) which include the static state WOB and 
the constraint loads resulting from drilling tool assembly contact with the wall and 
bottom of the wellbore, 
5 As the drilling tool assembly is incrementally "rotated**, constraint loads 

acting on the bit may change. For example, points of the drilling tool assembly in 
contact with the borehole surface prior to rotation may be moved along the surface 
of the wellbore resulting in friction forces at those points. Similarly, some points 
of the drilling tool assembly, which were nearly in contact with the borehole 

10 surface prior to the incremental rotation, may be brought into contact with the 
formation as a result of the incremental rotation, resulting in impact forces on the 
drilling tool assembly at those locations. As shown in Figure 7A-C, changes in the 
constraint loads resulting from the incremental rotation of the driUing tool 
assembly can be accounted for in the wall interaction update loop 245. 

15 In this example, once the system's response (i.e., new configuration) under 

the current loading conditions is obtained, the positions of the nodes in the new 
configuration are checked (at 244) in the wall constraint loop 245 to determine 
whether any nodal displacements fall outside of the bounds (Le,, violate constraint 
conditions) defined by the wellbore wall. If nodes are found to have moved 

20 outside of the wellbore wall, the impact and/or friction forces which would have 
occurred due to contact with the wellbore wall are approximated for those nodes 
(at 248) using the impact and/or friction models or parameters provided as input at 
208. Then the global load vector for the drilling tool assembly is updated (also 
shown at 208) to reflect the newly determined constraint loads. Constraint loads 

25 to be calculated may be determined to result from impact if, prior to the 
incremental rotation, the node was not in contact with the wellbore wall. 
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Similarly, the constraint load can be determined to result from frictional drag if the 
node now in contact with the wellbore wall was also in contact with the wall prior 
to the incremental rotation. Once the new constraint loads are determined and the 
global load vector is updated, at 248, the drilling tool assembly response is 

5 recalculated (at 244) for the same incremental rotation under the newly updated 
load vector (as indicated by loop 245). The nodal displacements are then 
rechecked (at 246) and the wall interaction update loop 245 is repeated until a 
dynamic response within the wellbore constraints is obtained. 

Once a dynamic response conforming to the borehole wall constraints is 

10 deteimined for the incremental rotation, the constraint loads on the drilling tool 
assembly due to interaction with the bottomhole during the incremental rotation 
are detemiined in the bit interaction loop 250. Those skilled in the art will 
appreciate that any method for modeling drill bit/earth formation interaction 
during drilling may be used to determine the forces acting on the drill bit during 

1 5 the incremental rotation of the drilling tool assembly. An example of one method 
is illustrated in the bit interaction loop 250 in Figure 7A-C. 

In the bit interaction loop 250, the mechanics analysis model of the drill bit 
is subjected to the incremental rotation angle calculated for the lowest node of the 
drilling tool assembly, and is then moved laterally and vertically to the new 

20 position obtained from the same calculation, as shown at 249. As previously 
noted, the drill bit in this example is a roller cone drill bit. Thus, in this example, 
once the bit rotation and new bit position are determined, interaction between each 
cone and the formation is determined. For a first cone, an incremental cone 
rotation angle is calculated at 252 based on a calculated incremental cone rotation 

25 speed and used to determine the movement of the cone during the incremental 
rotation. It should be understood that the incremental cone rotation speed can be 
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detennined from all the forces acting on the cutting elements of the cone and 
Newton's second law of motion. Altematively, it may be approximated from the 
rotation speed of the bit and the effective radius of the "drive row" of the cone. 
The effective radius is generally related to the lateral extent of the cutting elements 
5 that extend the farthest from the axis of rotation of the cone. Thus, the rotation 
speed of the cone can be defined or calculated based on the calculated bit 
rotational speed and the defined geometry of the cone provided as input (e.g., the 
cone diameter profile, cone axial offset, etc.) 

Then, for the first cone, interaction between each cutting element and the 

10 earth formation is determined in the cutting element/formation interaction loop 
256. In this interaction loop 256, the new position of a cutting element, for 
example, cutting element j on row ky is calculated 258 based on the incremental 
cone rotation and bit rotation and translation. Then, the location of cutting 
element j\k relative to the bottomhole and wall of the wellbore is evaluated, at 259, 

15 to determine whether cutting element interference (or contact) with the formation 
occurred during the incremental rotation of the bit If it is determined that contact 
did not occur, tiien the next cutting element is analyzed and the interaction 
evaluation is repeated for the next cutting element. If contact is determined to 
have occurred, then a depth of penetration, interference projection area, and 

20 scraping distance of the cutting element in the formation are determined, at 262, 
based on the next movement of the cutting element during the incremental 
rotation. The depth of penetration is the distance from the earth formation surface 
a cutting element penetrates into the earth fomiation. Depth of penetration can 
range from zero (no penetration) to the fall height of the cutting element (fiill 

25 penetration). Interference projection area is the fractional amount of the cutting 
element surface area, corresponding to the depth of penetration, which actually 
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contacts the earth formation. A fractional amount of contact usually occurs due to 
craters in the formation formed from previous contact with cutting elements. 
Scraping distance takes into account the movement of the cutting element in the 
formation during the incremental rotation. Once the depth of penetration, 

5 interference projection area, and scraping distance are determined for cutting 
element jjk these parameters are used in conjunction with - the cutting 
element/formation interaction data to detennine the resulting forces (constraint 
forces) exerted on the cutting element by the earth formation (also indicated at 
262). For example, force may be determined using the relationship disclosed in 

10 U.S. Patent Application No. 09/524,088, noted above and incorporated herein by 
reference. 

Once the cutting element/formation interaction variables (area, depth, force, 
etc.) are determined for cutting element jX the geometry of the bottom surface of 
the wellbore can be temporarily updated, at 264, to reflect the removal of 

15 formation by cutting element jjc during tiie incremental rotation of the drill bit. 
The actual size of the crater resulting from cutting element contact with the 
formation can be determined from the cutting element/earth formation interaction 
data based on the bottomhole surface geometry, and the forces exerted by the 
cutting element. One such procedure is described in U.S. Patent Application No. 

20 09/524,088, noted above. 

After the bottomhole geometry is temporarily updated, insert wear and 
strength can also be analyzed, as shown at 270, based on wear models and 
calculated loads on tiie cutting elements to deteraiine wear on the cutting elements 
resulting from contact with the formation and tiie resulting reduction in cutting 

25 element sfa-engtii. Then, the cutting element/formation interaction loop 260 
calculations are repeated for tiie next cutting element (/=y+l) of row k until cutting 



28 



element/formation interaction for each cutting element of tiie row is determined. 

Once the forces on each cutting element of a row are determined, the total 
forces on that row are calculated (at 268) as a sum of all the forces on the cutting 
elements of that row. Then, the cutting element/earth formation interaction 
5 calculations are repeated for the next row on the cone (^Ar+l) (in tiie row 
interaction loop 269) until the forces on each of the cutting elements on each of 
the rows on that cone are obtained. Once interaction of all of the cutting elements 
on a cone is determined, cone shell interaction with the formation is determined by 
checking node displacements at the cone surface, at 270, to determine if any of the 

10 nodes are out of bounds with respect to (or make contact with) the wellbore wall 
or bottomhole surface. If cone shell contact is determined to have occurred for the 
cone during tfie incremental rotation, the contact area and depth of penetration of 
the cone shell are determined (at 272) and used to determine interaction forces on 
the cone shell resulting from the contact. 

1 5 Once forces resulting from cone shell contact with the formation during the 

incremental rotation are determined, or it is determined that no shell contact has 
occurred, the total interaction forces on the cone during the incremental rotation 
can be calculated by summing all of the row forces and any cone shell forces on 
the cone, at 274. The total forces acting on the cone during the incremental 

20 rotation may then be used to calculate the incremental cone rotation speed ^, at 
276. Cone interaction calculations are then repeated for each cone (/=/+!) until 
the forces, rotation speed, etc. on each of the cones of the bit due to interaction 
with the formation are determined. 

Once the mteraction forces on each cone are determined, the total axial 

25 force on the bit (dynamic WOB) during the incremental rotation of the drilling 
tool assembly is calculated 278, from the cone forces. The newly calculated bit 



29 



interaction forces are then used to update the global load vector (at 279), and the 
response of the drilling tool assembly is recalculated (at 280) under the updated 
loading condition. The newly calculated response is then compared to the 
previous response (at 282) to determine if the responses are substantially similar. 

5 If the responses are determined to be substantially similar, then the newly 
calculated response is considered to have converged to a correct solution. 
However, if the responses are not determined to be substantially similar, then the 
bit interaction forces are recalculated based on the latest response at 284 and the 
global load vector is again updated (as indicated at 284). Then, a new response is 

1 0 calculated by repeating the entire response calculation (including the wellbore wall 
constraint update and drill bit interaction force update) until consecutive responses 
are obtained which are determined to be substantially similar (indicated by loop 
285), thereby indicating convergence to the solution for dynamic response to the 
incremental rotation. 

1 5 Once the dynamic response of the drilling tool assembly to an incremental 

rotation is obtained from the response force update loop 285, the bottomhole 
surface geometry is then permanently updated (at 286) to reflect the removal of 
formation corresponding to the solution. At this point, output information desired 
from the incremental simulation step can be provided as output or stored For 

20 example, the new position of the drilling tool assembly, the dynamic WOB, cone 
forces, cutting element forces, inipact forces, friction forces, may be provided as 
output information or stored. 

This dynamic response simulation loop 240 as described above is then 
repeated for successive incremental rotations of the bit until an end condition of 

25 the simulation (checked at 290) is satisfied. For example, using the total number 
of bit revolutions to be simulated as the termination command, the incremental 
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rotation of the drilling tool assembly and subsequent iterative calculations of flie 
dynamic response simulation loop 240 will be repeated until the selected total 
number of revolutions to be simulated is reached. Repeating the dynamic 
response simulation loop 240 as described above will result in simulating the 
5 performance of an entire drilling tool assembly drilling earth formations with 
continuous updates of the bottomhole pattern as drilled, thereby simulating the 
drilling of the drilling tool assembly in the selected earth formation. Upon 
completion of a selected number of operations of the dynamic response simulation 
loop, results of the simulation may be used to generate output information at 294 

10 characterizing the performance of the drilling tool assembly drilling the selected 
earth formation under the selected drilling conditions, as shown in Figure 7A-C. It 
should be understood that the simulation can be stopped using any other suitable 
tCTmination indicator, such as a selected wellbore depth desired to be drilled, 
indicated divergence of a solution, etc. 

15 As noted above, output information from a dynamic simulation of a drilling 

tool assembly drilling an earth formation may include, for example, the drilling 
tool assembly configuration (or response) obtained for each time increment, and 
corresponding bit forces, cone forces, cutting element forces, impact forces, 
friction forces, dynamic WOB, resulting bottomhole geometry, etc. This output 

20 information may be presented in the form of a visual representation (indicated at 
294), such as a visual representation of the borehole being drilled tiirough the earth 
formation with continuous updated bottomhole geometries and the dynamic 
response of the drilling tool assembly to drilling presented on a computer screen. 
Alternatively, tiie visual representation may include graphs of parameters provided 

25 as input and/or calculated during the simulation. For example, a time history of 
the dynamic WOB or the wear of cutting elements during drilling may be 
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presented as a graphic display on a computer screen. It should be understood that 
the invention is not limited to any particular type of display. Further, the means 
used for visually displaying aspects of simulated drilling is a matter of 
convenience for the system designer, and is not intended to limit the invention. 

5 One example of output data converted to a visual representation is illustrated in 
Figure 12, wherein the rotation of the drilling tool assembly and corresponding 
drilling of the formation is graphically illustrated as a visual display of drilling and 
desired parameters calculated during drilling can be numerically displayed. 

The example described above represents only one embodiment of the 

10 invention. Those skilled in the art will appreciate that other embodiments can be 
devised which do not depart from the scope of the invention as disclosed herein. 
For example^ an alternative method can be used to account for changes in 
constraint forces during incremental rotation. For example, instead of using a 
fmite element method, a finite difference method or a weighted residual method 

15 can be used to model the drilling tool assembly. Similarly, other methods may be 
used to predict the forces exerted on the bit as a result of bit/cutting element 
interaction with the bottomhole surface. For example, in one case, a method for 
interpolating between calculated values of constraint forces may be used to predict 
the constraint forces on the drilling tool assembly or a different method of 

20 predicting the value of the constraint forces resulting firom impact or fiictional 
contact may be used. Further, a modified version of the method described above 
for predicting forces resulting fi-om cutting element interaction with the 
bottomhole surface may be used. These methods can be analytical, numerical 
(such as finite element method), or experimental. Altematively, methods such as 

25 disclosed in SPE Paper No. 29922 noted above or PCT Patent Application Nos. 
WO 00/12859 and WO 00/12860 may be used to model roller cone drill bit 
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interaction with the bottomhole surface, or methods such as disclosed in SPE 
papers no. 15617 and no. 15618 noted above may be used to model fixed-cutter bit 
interaction with the bottomhole surface if a fixed-cutter bit is used. 



5 Method for Designing a Drilling Tool Assembly 

In another aspect, the invention provides a method for designing a drilling 
tool assembly for drilling earth formations. For example, the method may include 
simulating a dynamic response of a drilling tool assembly, adjusting the value of at 
least one drilling tool assembly design parameter, repeating the simulating, and 

10 repeating the adjusting and the simulating until a value of at least one drilling 
performance parameter is determined to be an optimal value; 

Methods in accordance with this aspect of the invention may be used to 
analyze relationships between drilling tool assembly design parameters and 
drilling performance of a drilling tool assembly. This method also may be used to 

15 design a drilling tool assembly having enhanced drilling characteristics. Further, 
the method may be used to analyze the effect of changes in a drilling tool 
configuration on drilling performance. Additionally, the method may enable a 
drilling tool assembly designer or operator to determine an optimal value of a 
drilling tool assembly design parameter for drilling at a particular depth or in a 

20 particular formation. 

Examples of drilling tool assembly design parameters include the type and 
number of components included in the drilling tool assembly; the length, ID, OD, 
weight, and material properties of each component; and the type, size, weight, 
configuration, and material properties of the drill bit; and the type, size, number, 

25 location, orientation, and material properties of the cutting elements on tihe bit. 
Material properties in designing a drilling tool assembly may include, for example, 
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the strength, elasticity, and density of the material. It should be understood that 
drilling tool assembly design parameters may include any other configuration or 
material parameter of the drilling tool assembly without departing from the spirit 
of the invention. 

5 Examples of drilling performance parameters include rate of penetration 

(ROP), rotary torque required to turn the drilling tool assembly, rotary speed at 
which the drilling tool assembly is turned, drilling tool assembly vibrations 
induced during drilling (e.g., lateral and axial vibrations), weight on bit (WOB), 
and forces acting on the bit, cones, and cutting elements. Drilling performance 

10 parameters may also include the inclination angle and azimuth direction of the 
borehole being drilled. One skilled in the art will appreciate that other drilling 
performance parameters exist and may be considered as detennined by the drilling 
tool assembly designer without departing from the spirit of the invention. 

In one appUcation of this aspect of the invention, illustrated in Figure 8, the 

15 method comprises defining, selecting or otherwise providing initial input 
parameters at 300 (including drilling tool assembly design parameters). The 
method further comprises simulating the dynamic response of the drilling tool 
assembly at 310, adjustmg at least one drilling tool assembly design parameter at 
320, and repeating the simulating of the drilling tool assembly 330. The method 

20 also comprises evaluating the change in value of at least one drilling performance 
parameter 340, and based on that evaluation, repeating the adjusting, the 
simulating, and the evaluating until at least one drilling performance parameter is 
optimized. 

As shown in the more detailed example of Figure 9, the initial parameters 
25 400 may include initial drilling tool assembly parameters 402, initial drilling 
environment parameters 404, drilling operating parameters 406, and drilling tool 
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assembly/drilling environment interaction parameters and/or models 408. These 
parameters may be substantially the same as the input parameters described above 
for the previous aspect. 

In this example, simulating 411 comprises constructing a mechanics 
5 analysis model of the drilling tool assembly (at 412) based on the drilling tool 
assembly parameters 402, determining system constraints at 414 using the drilling 
environment parameters 404, and then using the mechanics analysis model along 
with the system constraints to solve for the initial static state of the drilling tool 
assembly in flie drilling environment (at 416). Simulating 411 fiirther comprises 

10 using the mechanics analysis model along with the constraints and drilling 
operation parameters 406 to incrementally solve for the response of the drilling 
tool assembly to rotational input from a rotary table (at 418) and/or downhole 
motor, if used. In solving for the dynamic response, the response is obtained for 
successive incremental rotations until an end condition signaling the end of the 

1 5 simulation is detected. 

Incrementally solving for the response may also include determining, from 
drilling tool assembly/environment interaction information, loads on the drilling 
tool assembly during the incremental rotation resulting from changes in interaction 
between the drilling tool assembly and the drilling environment during the 

20 incremental rotation, and then recalculating the response of the drilling tool 
assembly under the new constraint loads* Incrementally solving may further 
include repeating, if necessary, the deteraiining loads and the recalculating of the 
response until a solution convergence criterion is satisfied. 

Examples for constracting a mechanics analysis model, detemiining initial 

25 system constraints, determining the initial static state, and incrementally solving 
for the dynamic response of the drilling tool assembly are described in detail for 
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the previous aspect of the invention. 

In the present example shown in Figure 9, adjusting at least one drilling 
tool assembly design parameter 426 comprises changing a value of at least one 
drilling tool assembly design parameter after each simulation by data iiq)ut j&om a 
5 file, data input from an operator, or based on calculated adjustment factors in a 
simulation program, for example. 

Drilling tool assembly design parameters may include any of the drilling 
tool assembly parameters noted above. Thus in one example, a design parameter, 
such as the length of a drill collar, can be repeatedly adjusted and simulated to 

10 determine the effects of BHA weight and length on a drilling performance 
parameter ROP). Similarly, the inner diameter or outer diameter of a drilling 
collar may be repeatedly adjusted and a corresponding change response obtained. 
Similarly, a stabilizer or other component can be added to the BHA or deleted 
from the BHA and a corresponding change in response obtained. Further, a bit 

15 design parameter may be repeatedly adjusted and corresponding dynamic 
responses obtained to determine the effect of changing one or more drill bit design 
parameters, such as cone profile, insert shape and size, number of rows offsets (for 
roller cone bits) on the drilling performance of the drilling tool assembly. 

In the example of Figure 9, repeating the simulating 411 for the "adjusted" 

20 drilling tool assembly comprises constructing a new (or adjusted) mechanics 
analysis model (at 412) for the adjusted drilling tool assembly, determining new 
system constraints (at 414), and then using the adjusted mechanics analysis model 
along with Ae corresponding system constramts to solve for the initial static state 
(at 416) of the of the adjusted drilling tool assembly in the drilling environment. 

25 Repeating the simulating 411 further comprises using the mechanics analysis 
model, initial conditions, and constraints to incrementally solve for the response of 
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the adjusted drilling tool assembly to simulated rotational input from a rotary table 
and/or a downhole motor, if used. 

Once the response of the previous assembly design and the response of the 
current assembly design are obtained, the effect of the change in value of at least 
5 one design parameter on at least one drilling performance parameter can be 
evaluated (at 422). For example, during each simulation, values of desired drilling 
performance parameters (WOB, ROP, impact loads, etc) can be calculated and 
stored. Then, these values or other factors related to the drilling response (such as 
vibration factors), can be analyzed to determine the effect of adjusting the drilling 

10 tool assembly design parameter on the value of the at least one drilling 
performance parameter. 

Once an evaluation of at least one drilling parameter is made, based on that 
evaluation the adjusting and the simulating may be repeated until it is determined 
that the at least one drilling performance parameter is optimized or an end 

15 condition for optimization has been reached (at 424). A drilling performance 
parameter may be determined to be at an optimal value when a maximum rate of 
penetration, a minimum rotary torque for a given rotation speed, and/or most even 
weight on bit is determine for a set of adjustment variables. Other drilling 
performance parameters, such as minimized lateral impact force or 

20 optimized/balanced forces on different cones for roller cone bit applications can 
also be used. A simplified example of repeating the adjusting and the simulating 
based on evaluation of consecutive responses is as follows. 

Assume that the BHA weight is the drilling tool assembly design parameter 
to be adjusted (for example, by changing the length, equivalent ID, OD, adding or 

25 deleting components), and ROP is the drilling performance parameter to be 
optimized. Therefore, after obtaining a first response for a given drilling tool 
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assembly configuration, the weight of the BHA can be increased and a second 
response can be obtained for the adjusted drilling tool assembly. The weight of 
the BHA can be increased, for example, by changing the ID for a given OD of a 
collar in the BHA (will ultimately affect the system mass matrix). Alternatively, 

5 the weight of the BHA can be increased by increasing the length, OD, or by 
adding a new collar to the BHA (will ultimately affect the system stiffiiess matrix). 
In either case, changes to the drilling tool assembly will effect the mechanics 
analysis model for the system and the resulting initial conditions. Therefore, the 
mechanics analysis model and initial conditions will have to be re-determined for 

10 the new configuration before a solution for the second response can be obtained. 
Once the second response is obtained, the two responses (one for the old 
configuration, one for the new configuration) can be compared to determine which 
configuration (BHA weight) resulted in the most favorable (or greater) ROP. If 
the second configuration is found to result in a greater ROP, then the weight of the 

1 5 BHA may be further increased, and a (third) response for the newer configuration) 
may be obtained and compared to the second. Alternatively, if the increase in the 
weight of the BHA is found to result in a decrease in the ROP, then the drilling 
tool assembly design may be readjusted to decrease the BHA weight to a value 
lower than that set for the first drilling tool assembly configuration and a (third) 

20 response may be obtained and compared to the first. This adjustment, 
recalculation, evaluation may be repeated until it is detemiined that an optimal or 
desired value of at least one drilling performance parameter, such as ROP in this 
case, is obtained. 

Advantageously, embodiments of the invention may be used to analyze the 
25 relationship between drilling tool assembly design parameters and drilling 
performance in a selected drilling environment. Additionally, embodiments of the 
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invention may be used to design a drilling tool assembly having optimal drilling 
performance for a given set of drilling conditions. Those skilled in the art will 
appreciate that other embodiments of the invention exist which do not depart from 
the spirit of this aspect of the invention. 

5 

Method for Qptimiying Drilling Operating Parameters for a Selected- or Particular 
prillinp Tool Assembly 

In another aspect, the invention provides a method for determining optimal 
drilling operating parameters for a selected drilling tool assembly. In one 
1 0 embodiment, this method includes simulating a dynamic response of a drilling tool 
assembly, adjusting the value of at least one drilling operating parameters, 
repeating the simulating, and repeating the adjusting and the simulating until a 
value of at least one drilling performance parameter is determined to be an optimal 
value. 

1 5 The method in accordance with this aspect of the invention may be used to 

analyze relationships between drilling operating parameters and the drilling 
performance of a selected drilling tool assembly. The method also may be used to 
improve the drilling performance of a selected drilling tool assembly. Further, the 
method may be used to analyze the effect of changes in drilling operating 

20 parameters on the drilling performance of the selected drilling tool assembly. 
Additionally, the method in accordance with this aspect of the invention may 
enable the drilling tool assembly designer or operator to determine optimal drilling 
operating parameters for a selected drilling tool assembly drilling a particular 
depth or in a particular formation. 

25 As previously explained, drilling operating parameters include, for 

example, rotational speed at which the drilling tool assembly is turned, or rotary 
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torque applied to turn the drilling tool assembly, hook load (which is one of the 
major factors to influence WOB), drilling fluid flow rate, and material properties 
of the drilling fluid (e.g., viscosity, density, etc). It should be understood that 
drilling parameters may include any drilling environment or drilling operating 
5 parameters which may affect the drilling performance of a drilling tool assembly 
without departing from the spirit of the invention- 
Drilling performance parameters that may be considered in optimizing the 
design of a driUing tool assembly may include, for example, the ROP, rotary 
torque required to turn the drilling tool assembly, rotary speed at which the drilling 

10 tool assembly is turned, drilling tool assembly vibrations (in terms of velocities, 
accelerations, etc.), WOB, lateral force, moments, etc. on the bit, lateral and axial 
forces, moments, etc. on the cones, and lateral and axial forces on the cutting 
elements. It should be understood that during simulation velocity and 
displacement are calculated for each node point and can be used to calculate 

15 force/acceleration as an indicator of drilling tool assembly vibrations. One skilled 
in the art will appreciate that other parameters which can be used to evaluate 
drilling performance exist and may be used as determined by the drilling tool 
assembly designer without departing from the spirit of the invention. 

Figure 10 shows a flow chart for one example of a method for determining 

20 at least one optimal drilling operating parameter for a selected drilling tool 
assembly. In this example, the method comprises defining, selecting or otherwise 
providing initial input parameters at 500 (including drilling tool assembly design 
parameters and drilling operating parameter) which describe various aspects of the 
initial system. The method further comprises simulating the dynamic response of 

25 a drilling tool assembly at 510, adjusting at least one drilling operating parameter 
at 520, and repeating the simulating of the drilling tool assembly at 530. The 
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method also comprises evaluating the change in value of at least one drilling 
performance parameter 540, and based on that evaluation, repeating the adjusting 
520, the simulating 530, and the evaluating 540 until at least one drilling 
performance parameter is optimized. 

5 Another example of such a method is shown in Figure 11. In this example, 

the initial parameters 600 include initial drilling tool assembly parameters 602, 
initial drilling environment parameters 604, initial drilling operating parameters 
606, and drilling tool assembly/drilling environment interaction parameters and/or 
models 608. These parameters may be substantially the same as those described 

1 0 for the first aspect of the invention discussed above. 

In this example, once the input parameters 600 are provided, the input 
parameters 600 are used to construct a mechanics analysis model (at 612) of the 
drilling tool assembly and used to determine system constraints (at 614) (wellbore 
wall and bottom surface constraints). Then, the mechanics analysis model and 

15 system constraints are used to determine the initial conditions (at 616) oii the 
drilling tool assembly inserted in the wellbore. Examples for constructing a 
mechanics analysis model of a drilling tool assembly and determining initial 
constraints and initial conditions are described in detail above for the first aspect 
of the invention. 

20 In the example shown in Figure 11, simulatmg the dynamic response 611 

comprises using the mechanics analysis model along with the initial constraints 
and initial conditions to incrementally solve for the dynamic response of the 
drilling tool assembly to simulated rotational input firom a rotary table (at 618) 
and/or downhole motor. The dynamic response to successive incremental 

25 rotations is incrementally obtained until an end condition signaling the end of the 
simulation is detected. 
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Incrementally solving for the response may include iteratively determining, 
from drilling tool assembly/environment interaction data or models, new drilling 
environment interaction forces on the drilling tool assembly resulting from 
changes in interaction between the drilling tool assembly and the drilling 
environment during the incremental rotation, and then recalculating the response 
of the drilling tool assembly to the incremental rotation under the newly calculated 
constraint loads. Incrementally solving may further include repeatmg, if 
necessary, the determining and the recalculating until a constraint load 
convergence criterion is satisfied. An example of incrementally solving for the 
response as described here is presented in detail for the first aspect of the 
invention. 

At least one drilling operating parameter may be adjusted (at 620) as 
discussed above for the previous aspect of the invention, such as by reading in a 
new value from a data file, data input from an op^ator, or calculating adjustment 
valuejs based on evaluation of responses corresponding to previous values, for 
example. Similarly, drilling performance parameter(s) adjusted may be any 
parameter effecting the operation of drilling without departing from the spirit of 
the invention. In some cases, adjusted drilling parameters may be limited to only 
particular parameters. For example, the drilling tool assembjy designer/operator 
may concentrate only on the effect of the rotary speed and hook load (or WOB) on 
drilling performance, in which case only parameters effecting the rotary speed or 
hook load (or WOB) may be adjustable. 

In the example shown in Figure 11, rqpeating the simulating 618 comprises 
at least recalculating the response of the drilling tool assembly to the adjusted 
drilling operating conditions. However, if an adjustment is made to a drilling 
operating parameter that affects the drilling environment, such as the viscosity or 
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density of drilling fluid, repeating the simulation may comprise first determining a 
new system global damping matrix and global load vectors and then using the 
newly updated mechanics analysis model to incrementally solve for the response 
of the drilling tool assembly to simulated rotation under the new drilling operating 
5 conditions. However, if the adjustinent made to a drilling operating parameters 
does not affect the drilling environment, which may typically be the case (e.g., 
rotation speed of the rotary table), repeating the simulation may only comprise 
solving for the dynamic response of the drilling tool assembly to the adjusted 
operating conditions and the same initial conditions (the static equilibrium state) 
10 by using the mechanics analysis model. 

Similar to the previous aspect, once a response for the previous adjusted 
operating parameters and a response for the current adjusted operating parameters 
are obtained, the effect the change in value of the drilling operating parameter on 
drilling performance can be evaluated (at 622). For example, during each 
1 5 simulation values of desired drilling performance parameters (WOB, ROP, impact 
loads, optimized force distribution on cutting elements, optimized/balanced for 
distribution on cones for roller cone bits, optimized force distribution on lades for 
PDC bits, etc.) can be calculated. Then, tiiese values or otiier factors related to the 
response (such as vibration parameters) can be analyzed to determine the effect of 
20 adjusting the drilling operating parameter on the value of at least one drilling 
performance parameter. 

Optimization criteria may include optimizing the force distribution on 
cutting elements, maximizing the rate of penetration (ROP), minimizing the WOB 
required to obtain a given ROP, minimizing lateral fanpact force, etc. In addition, 
25 for roller cone drill bits, optimization criteria may also include optimizing or 
balancing force distribution on cones. For fixed-cutter bits, such as PDC bits, 
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optimization criteria may also include optimizing force distribution on the blades 
or among the blades. 

Once an evaluation of the least one drilling operating parameter is made, 
based on that evaluation the adjusting and the simulating may be repeated until it 
is determined that at least one drilling performance parameter is optimized, or 
until an end condition for optimization is reached. As noted for the previous 
aspect, a drilling performance parameter may be determined to be at an optimal 
value when, for example, a maximum rate of penetration, a minimum rotary 
torque for a given rotation speed, and/or most even weight on bit is detennine for a 
set of adjustment variables. Additionally, an end condition for optimization may 
include detennining when a change in the operation value no long results in an 
improvement in the drilling performance of the drilling tool assembly. A 
simplified example of repeating the adjusting, the simulating, and the evaluating 
until a drilling performance parameter is optimized is as follows. 

For example, if after obtaining a first response, the hook load is decreased 
(which ultimately increases the WOB), and then a second response is obtained for 
the decreased hook load, the ROP of the two responses can be compared. If the 
second response is found to have a greater ROP than the first (/.e., decreased hook 
load is shown to increase ROP), the hook load may be fiirther decrease and a third 
response may be obtained and compared to the second. This adjustment, 
resimulation, evaluation may be repeated until the point at which decrease in hook 
load provides maximum ROP is obtained. Alternatively, if the decrease in hook 
load is found to result in an decrease in the ROP, then the hook load may be 
increased to value higher than the value of the hook load for the first simulation, 
and a third response may be obtained and compared with the first (having the more 
favorable ROP). This adjustment, resimulation, evaluation may be repeated until 
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it is detennined that further increase in hook load provides no further benefit in the 
ROP. 

Advantageously, embodiments of the invention may be used to analyze the 
relationship between drilling parameters and drilling performance for a select 
5 drilling tool assembly drilling a particular earth formation. Additionally, 
embodiments of the invention may be used to optimize the drilling performance of 
a given drilling tool assembly. Those skilled in the art will appreciate that other 
embodiments of the invention exist which do not depart from the spirit of this 
aspect of the invention. 

10 Further, it should be understood that regardless of the complexity of a 

drilling tool assembly or the trajectory of the wellbore in which it is to be 
constrained, the invention provides reliable methods that can be used for 
predicting the dynamic response of the drilling tool assembly drilling an earth 
fonmation. The invention also facilitates designing a drilling tool assembly having 

15 enhanced drilling performance, and helps deteimine optimal drilling operating 
parameters for improving the drilling performance of a selected drilling tool 
assembly. 

While the invention has been described with respect to a limited number of 
embodiments and examples, those skilled in the art will appreciate that other 
20 embodiments can be devised which do not depart from the scope of the invention 
as disclosed herein. Accordingly, the scope of the invention should be limited 
only by the attached claims. 
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CLAIMS 



1 1 . A method for optimizing a drilling tool assembly design, comprising: 

2 simulating a dynamic response of the drilling tool assembly; 

3 adjusting a value of at least one drilling tool assembly design parameter; 

4 repeating the simulating; 

5 repeating the adjusting and the simulating until at least one drilling 

6 performance parameter is determined to be at an optimal value. 

1 2. The method of claim 1, wherein the simulating comprises, 

2 solving for the dynamic response of the drilling tool assembly to an 

3 incremental rotation using a mechanics analysis model, and 

4 repeating said solving for a select number of successive incremental 

5 rotations. 

1 3 . The method of claim 2, wherein said solving comprises, 

2 constructing the mechanics analysis model of the drilling tool assembly 

3 using selected drilling tool assembly design parameters, 

4 determining wellbore constraints from wellbore trajectory parameters, a 

5 specified bottom hole geometry, and a specified hook load, 
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6 determining loads on the drilling tool assembly for a position of the drilling 

7 tool assembly in the wellbore using at least the mechanics analysis model and the 

8 wellbore constraints, and 

9 calculating the dynamic response of the drilling tool assembly under the 
1 0 loads to the incremental rotation using the mechanics analysis model. 

1 4. The method of claim 3, wherein said solving further comprises, 

2 redetermining the loads on the drilling tool assembly based on the 

3 calculated dynamic response to the incremental rotation, 

4 repeating the calculating the dynamic response of the drilling tool assembly 

5 under the loads to the incremental rotation, and 

6 repeating the redetermining and the calculating until convergence of the 

7 dynamic response is determined. 

1 5. The method of claim 4, wherein redetermining the loads comprises, 

2 identifying, from the dynamic response, points along the drilling tool 

3 assembly which interact with the wellbore wall during the incremental rotation, 

4 determining, from drilling tool assembly/ environment interaction 

5 information, constraint forces at the points resulting from interaction with the 

6 wellbore wall, and 

7 updating the loads to include determined constraint forces. 

1 6. The method of claim 5, wherein redetermining the loads fiirther comprises, 

2 determining, from the dynamic response, drill bit parameters, and a drill bit 

3 model, cutting element interaction with a bottom of the wellbore during the 

4 incremental rotation, 
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5 determining, from drilling tool assembly/ environment interaction 

6 information, cutting element interaction, and the hook load, total forces on the bit 

7 resulting from the cutting element interaction with the bottom surface of the 

8 wellbore, and 

9 updating the loads to account for the newly calculated total forces on the 
10 bit. 

1 7. The method ofclaim 3, wherein the determining loads comprises, 

2 determining constraint forces required to displace the drilling tool assembly 

3 from an unconstrained state to a state wherein a centerline of the drilling tool 

4 assembly substantially aligns with a centerline of a wellbore trajectory, 

5 calculating the steady state position of the drilling tool assembly under the 

6 determined constraint forces, 

7 redetermining constraint forces required to constrain the steady state 

8 position of the drilling tool assembly within the wellbore, and 

9 repeating the calculating of Uie steady state position and the redetermining 
10 of the constraint forces until a position convergence criterion is satisfied. 

1 8. The method of claim 1, wherein the drilling performance parameter is 

2 determined to be at an optimal value when at least one of a maximum rate of 

3 penetration, a minimum rotary torque to maintain rotation speed, and a most even 

4 weight on bit is determined to occur. 

1 9. The method of claim 1, wherein the at least one drilling tool assembly 

2 design parameter is selected from the group of drill string design parameters, 

3 bottomhole assembly, design parameters, and drill bit design parameters. 
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1 10. The method of claim 9, wherein the drill string design parameters comprise 

2 at least one of a length, an inner diameter, an outer diameter, a density, a strength, 

3 and an elasticity for at least one component in a drill string, wherein the drill string 

4 conq)rises at least one joint of drill pipe. 

1 11. The method of claim 9, wherein the bottomhole assembly design 

2 parameters comprise at least one selected from the group of a length, an inner 

3 diameter, an outer diameter, a weight, a strength, and an elasticity for at least one 

4 of a plurality of components in a bottomhole assembly, adding at least one 

5 component to the bottomhole assembly, and deleting at least one component from 

6 the bottomhole assembly, wherein components in the bottomhole assembly 

7 comprise at least one of a drill collar, stabilizer, bent housing, measurement-while- 

8 drilling tool, logging-while-drilling tool, and downhole motor. 

1 12. The method of claim 9, wherein the drill bit design parameters comprise at 

2 least one of a drill bit type, drill bit diameter, cutting element count, cutting 

3 element geometric shape, cutting element height, cutting element location, and 

4 cutting element spacing. 

1 13. The method of claim 12, wherein the drill bit type is a roller cone drill bit, 

2 and the drill bit design parameters further comprise at least one of a number of 

3 cones, a cone profile, a number of cutting element rows on each cone, a number of 

4 cutting elements on each row, a cutting element orientation, a cutting element 

5 pitch, a cone axis offeet, and a journal angle. 
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14. The method of claim 1, wherein the at least one drilling performance 
parameter is selected jfrom the group of rate of penetration, rotary torque, rotary 
speed, weight on bit, lateral force on bit, ratio of forces on cones, ratio of forces 
between cones, distribution of forces on cutting elements, volume of formation 
cut, and wear on cutting elements. 

1 5. A method for determining at least one optimal drilling operating parameter 
for a drilling tool assembly, comprising: 

simulating a dynamic response of the drilling tool assembly; 
adjusting a value of at least one drilling operating parameter, 
repeating the simulating; 

repeating the adjusting and the simulating until at least one drilling 
performance parameter is determined to be at an optimal value. 

1 6. The method of claim 1 5, wherein the simulatmg comprises, 

solving for the dynamic response of the drilling tool assembly to an 
incremoital rotation using a mechanics analysis model, and 

repeating said solving for a select number of successive incremental 
rotations. 

1 7. The method of claim 1 6, wherein said solving comprises 
determining wellbore constraints on the drilling tool assembly, 
determining loads on the drilling tool assembly resulting from wellbore 

constraints, and 

calculating the dynamic response of the drilling tool assembly under the 
loads and drilling operating parameters to the incremental rotation. 
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1 1 8. The method of claim 1 7, wherein said solving further comprises, 

2 redetennining loads on the drilling tool assembly based on the dynamic 

3 response to the incremental rotation, 

4 repeating the calculating the dynamic response of the drilling tool assembly 

5 undCT the loads to the incremental rotation, and 

6 repeating the redetermining and calculating until convergence of the 

7 dynamic response is determined. 

1 19. The method of claim 1 8, wherem redetennining the loads comprises, 

2 identifying, from the dynamic response, points along the drilling tool 

3 assembly which interact with the wellbore wall during the incremental rotation, 

4 determining, from drilling tool iassembly/ environment interaction 

5 information, constraint forces at the points resulting from interaction with the 

6 wellbore wall, and 

7 updating the loads to include determined constraint forces. 

1 20. The method of claim 19, wherein redetermining the loads further 

2 comprises, 

3 determining, from the dynamic response, drill bit parameters, and a drill bit 

4 model, cutting element interaction with a bottom of the wellbore during the 

5 incremental rotation, 

6 determining, from drilling tool assembly/ environment interaction 

7 infomiation, cutting element interaction, and the hook load, total forces on the bit 

8 resulting from the cutting element mteraction with the bottom surface of the 

9 wellbore, and 
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10 updating the loads to account for the newly calculated total forces on the 

11 bit. 

1 21. The method of claim 1 7, wherein the determining loads comprises, 

2 determining constraint forces required to displace the drilling tool assembly 

3 from an unconstrained state to a state wherein a centerline of the. drilling tool 

4 assembly substantially aligns with a centerline of a wellbore trajectory, 

5 calculating the steady state position of the drilling tool assembly under the 

6 determined constraint forces, 

7 redetermining constraint forces required to constrain the steady state 

8 position of the drilHng tool assembly within the wellbore, and 

9 repeating the calculating of the steady state position and the redetermining 
10 of the constraint forces until a position convergence criterion is satisfied. 



1 22. The method of claim 15, wherein the at least one drilling operating 

2 parameter is selected from the group of rotary speed, rotary torque, hook load, 

3 drilling fluid viscosity, and drilling fluid density. 

1 23. The method of claim 15, wherein the at least one drilling performance 

2 parameter is selected from the group of rate of penetration, rotary torque, rotary 

3 speed, weight on bit, lateral force on bit, ratio of forces on cones, distribution of 

4 forces on cutting elements, volume of formation cut, and wear on cutting 

5 elements. 

1 24. The method of claim 15, wherein the at least one drilling performance 

2 parameter is determined to be at an optimal value when at least one of a maximum 
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rate of penetration, a minimum rotary torque to maintain the rotation speed, and a 
most even weight on bit is determined to occur. 



1 25. A method for designing a driUing tool assembly, comprising: 

2 defining initial drilling tool assembly design parameters; 

3 simulating a dynamic response of the drilling tool assembly; ■ 

4 adjusting a value of at least one of the drilling tool assembly design 

5 parameters; 

6 repeating the simulating and the adjusting a selected number of times; 

7 evaluating the dynamic responses; and 

8 selecting, based on the evaluating, desired drilling tool assembly design 

9 parameters. 

1 26. A method for selecting drilling operating parameters for a drilling tool 

2 assembly, comprising: 

3 simulating a dynamic response of the drilling tool assembly; 

4 adjusting a value of at least one drilling operating parameter; 

5 repeating tihe simulating; 

6 repeating the adjusting and the simulating a selected number of times; 

7 evaluating the dynamic responses simulated; and 

8 selecting, based on the evaluating, drilhng operating parameter values. 

1 27. A method for generating a visual representation of a drilling tool assembly, 

2 drilling earth formation, the drilling tool assembly comprising at least a drill pipe 

3 and a drill bit, the method comprising: 
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solving for a dynamic response of the drilling tool assembly to an 
incremental rotation; 

determining, based on the dynamic response, parameters of craters removed 
from a bottomhole surface of the formation due to contact of the bit with the 
bottomhole surface during the incremental rotation; 

calculating a bottomhole geometry, wherein the craters are removed from 
the bottomhole surface; 

repeating said solving, determining, and calculating for a selected number 
of successive incremental rotations; and 

converting the dynamic response and the bottomhole geometry parameters 
into a visual representation. 

28. A method for generating a visual representation of a drilling tool assembly 
drilling an earth formation, comprising: 

selecting drilling tool assembly design parameters, comprising at least a 
length of drill pipe, a geometry of at least one cutting element on a drill bit, and a 
location of the at least cutting element; 

selecting drilling parameters, comprising at least a rotation speed of the 
drilling tool assembly and a wellbore bottomhole surface and; 

selecting an earth formation to be represented as drilled; 

calculating from said selected drilling tool assembly design parameters, 
said selected drilling parameters, and said earth formation, a dynamic response of 
the drilling tool assembly and a bottomhole geometry resulting from interaction 
between the at least one cutting element on the drill bit and the bottomhole 
surface; 
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14 incrementally rotating said drilling tool assembly, and repeating said 

15 calculating; and 

16 converting said drilling tool assembly parameters and said bottomhole 

17 geometry parameters into said visual representation. 
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